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Abstract 

The equations of in-medium gluodynamics are proposed. Their classical 
lowest order solution is explicitly shown for a color charge moving with con- 
stant speed. For nuclear permittivity larger than 1 it describes emission of 
Cherenkov gluons resembling results of classical electrodynamics. The choice 
of nuclear permittivity and Lorentz-invariance of the problem are discussed. 
Effects induced by the transversely and longitudinally moving (relative to the 
collision axis) partons at LHC energies are described. 

1 Introduction 

The collective effects observed in ultrarelativistic heavy-ion collisions at SPS 
and RHIC [H [2j [31 d] have supported the conjecture of quark-gluon plasma 
(QGP) formed in these processes. The properties and evolution of this medium 
are widely debated. At the simplest level it is assumed to consist of a set 
of current quarks and gluons. It happens however that their interaction is 
quite strong so that the notion of the strongly interacting quark-gluon plasma 
(sQGP) has been introduced. Moreover, this substance reminds an ideal liquid 
rather than a gas. Whether perturbative quantum chromodynamics (pQCD) 
is applicable to the description of the excitation modes of this matter is doubt- 
ful. Correspondingly, the popular theoretical approaches use either classical 
solutions of in-vacuum QCD equations at the initial stage or hydrodynamics 
at the final stage of its evolution. 

The collective excitation modes of the medium may however play a crucial 
role. One of the ways to gain more knowledge about the excitation modes is 
to consider the propagation of relativistic partons through this matter. Phe- 
nomenologically their impact would be described by the nuclear permittivity 
of the matter corresponding to its response to passing partons. Namely this 
approach is most successful for electrodynamical processes in matter. There- 
fore it is reasonable to modify QCD equations by taking into account collective 
properties of the quark-gluon medium. For the sake of simplicity we consider 
here the gluodynamics only. The generalization to quarks is straightforward. 



email: dremin@lpi.ru 



1 



The classical lowest order solution of these equations coincides with Abelian 
electrodynamical results up to a trivial color factor. One of the most spec- 
tacular of them is Cherenkov radiation and its properties. Now, Cherenkov 
gluons take place of Cherenkov photons [6] . Their emission in high energy 
hadronic collisions is described by the same formulae but with nuclear per- 
mittivity in place of the usual one. It should be properly defined. Actually, 
one considers them as quasiparticles, i.e. quanta of the medium excitations 
with properties determined by the permittivity. The interplay of medium 
properties and velocity of the particle is crucial for the radiation field. 

Another important problem of this approach is related to the notion of the 
rest system of the medium. The Lorentz invariance is lost if the permittivity 
is introduced^]. Therefore one has to choose the proper coordinate system 
where its definition is at work. While it is simple for macroscopic media in 
electrodynamics, one should consider partons moving in different directions 
with different energies in case of heavy-ion collisions. It has direct impact 
on properties of emitted particles. The fast evolution of the medium and its 
short lifetime differ it from common electrodynamical examples. 

All these problems are discussed in what follows. 

2 Equations of in-medium gluodynamics 

At the beginning let us remind the classical in- vacuum Yang-Mills equations 

= r, (i) 

F^ u = d^A" — d"A li — ig[A^, A u ], (2) 

where A» = A%T a ; A a (A° a = $ OJ A a ) are the gauge field (scalar and vector) 
potentials, the color matrices T a satisfy the relation [T a ,T&] = if a bcT c , = 
dfj, — ig[A/j,, •], J v {p, j) is a classical source current, H = c = 1 and the metric 
tensor is g^ v =diag(+, -,-,-). 

In the covariant gauge d^A^ = they are written as 

HA* = J fl + ig[A v ,d u A» + F v % (3) 

where □ is the d'Alembertian operator. It was shown [8] (and is confirmed 
in what follows) that in this gauge the classical gluon field is given by the 
solution of the corresponding Abelian problem. 

2 In principle, this deficiency is cured by the relativistic generalization of the notion of permit- 
tivity (e.g., see [7]). 
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The chrome-electric and chrome-magnetic fields are 



E" = F»°, (4) 

= i ( '"-'/•''-'. (5) 
or as functions of gauge potentials in vector notations 

dA 

E a = -grad$ a ~-gf+ gfabcAb^c, (6) 

B a = curlA a - -gf ab c[A b A c ]. (7) 
The equations of motion ([1]) in vector form are written as 

divE a - gf o6c A 6 E c = p a , (8) 
<9E 

curlB a - — S - 5 / afec ($ b E c + [A 6 B C ]) = j a . (9) 
at 

The Abelian equations of in-vacuum electrodynamics are obtained from 
Eq. ([3]) if the second term in its right-hand side is put equal to zero and 
color indices omitted. The medium is accounted if E is replaced by D = eE 
in F^ u , i.e. in Eq. Therefore the Eqs. ([!]), ([9]) in vector form are 

most suitable for their generalization to in-medium case. The equations of in- 
medium electrodynamics differ from in-vacuum ones by dielectric permittivity 
e 7^ 1 entering there as 

8 2 A 

AA-e^ = -j, (10) 

e(A* - e^) = -p. (11) 

The permittivity describes the matter response to the induced fields which 
is assumed to be linear and constant in Eqs. (|10p . (jlip . It is determined 
by the distribution of internal current sources in the medium. Then external 
currents are only left in the right-hand sides of these equations. 
Now, the Lorentz gauge condition is 

divA + e— = 0. (12) 

The Lorentz invariance is broken if e ^ 1 in front of the second terms in 
the left-hand sides. Then one has to deal within the coordinate system where 



3 e denotes the dielectric permittivity of the medium. The magnetic permittivity is put equal 
to 1 to simplify the formulae. 
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a substance is at rest. The values of e are determined just there. To cancel 
these requirements one must use Minkowski relations between D, E, B, H 
valid for a moving medium [7J. It leads to more complicated formulae, and 
we do not use them in this paper. 

The most important property of solutions of these equations is that while 
the in- vacuum (e = 1) equations do not admit any radiation processes, it 
happens for e / 1 that there are solutions of these equations with non-zero 
Poynting vector. 

Now we are ready to write down the equations of in-medium gluodynam- 
ics generalizing Eq. ([3]) in the same way as Eqs. (jlOH . (Ilip are derived in 
electrodynamics. We introduce the nuclear permittivity and denote it also by 
e since it will not lead to any confusion. After that one should replace E a in 
Eqs. ©, © by eE a and get: 

e(divE a - gf abc A b E c ) = p a , (13) 

curlB a - - gf abc (e<S> b B c + [A 6 B C ]) = j a . (14) 

The space-time dispersion of e is neglected here. 

In terms of potentials these equations are cast in the form: 

d 2 A 1 d 

A A a - e-jjj^ = -Ja ~ #/ & c (-(curl[A 6 , A c ] + [A 6 curlA c ]) + — (A fe $ c ) - 

dA 1 
e$&-7^- - e$ 6 grad$ c - -gf cmn [A b [A ]} + gef cmn $ b A m <S> n ),(l5) 



A *« " ^ = ~ + 5/a 6c (2A b gradd> c + A b ^ - e^ c ) + 

g 2 famnfnlb-^mM^b- (16) 

If the terms with explicitly shown coupling constant g are omitted, one gets 
the set of Abelian equations which differ from electrodynamical equations 
pop , (jlip by the color index a only. Their solutions are shown in the next 
section. The external current is ascribed to a parton fast moving relative to 
other partons "at rest". 

The potentials are linear in g because the classical current J M is linear 
also. Therefore omitted terms are of the order of g s and can be taken into 
account as a perturbation. It was done in j9j[T0j for in-vacuum gluodynamics. 
Here, the general procedure is the same. After getting explicit lowest order 
solution (see the next section) one exploits it together with the non-Abelian 
current conservation condition to find the current component proportional to 
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g 3 . Then with the help of Eqs. (|15p . (|16|) one finds the potentials up to the 
order g 3 . They can be represented as integrals convoluting the current with 
the corresponding in-medium Green function. The higher order corrections 
may be obtained in the same way. We postpone their consideration for further 
publications. 

The crucial distinction between Eq. © and Eqs. (fT5j) . (fTUj) is that there is 
no radiation (the field strength is zero in the forward light-cone and no gluons 
are produced) in the lowest order solution of Eq. ([3]) and it is admitted for Eqs. 
(fT5|) . (fT6j) because e takes into account the collective response (polarization) 
of the nuclear matter. We have assumed that no color indices are attached 
to e. It would correspond to the collective response of the color-neutral (on 
the average) medium if color exchange between the external current and 
medium excitations is numerous and averages to zero. The lack of knowledge 
about the collective excitations of the nuclear medium prevents more detailed 
studies. However it seems to be justified at least for Cherenkov effects. 

3 Cherenkov gluons as the classical lowest 
order solution of in-medium gluodynamics 

Cherenkov effects are especially suited for treating them by classical approach 
to Eqs. (fT5|) . (fT6|) . Their unique feature is independence of the coherence of 
subsequent emissions on the time interval between these processes. 

The problem of the coherence length for Cherenkov radiation was exten- 
sively studied [HJ [12]. It was shown that the w-component of the field of a 
current can be imitated by a set of oscillators with frequency uj situated along 
the trajectory. The waves from all oscillators add up in the direction given by 
the Cherenkov angle 6 independent on the length of the interval filled in by 
these oscillators. The phase disbalance Acj) between emissions with frequency 
co = k/y/e separated by the time interval At (or the length Az = vAt) is 
given by 

Acj) = uAt - kAzcosO = kAz(-^-= - cosO) (17) 

v^Jt 

up to terms which vanish for large distances between oscillating sources and 
the detector. For Cherenkov effects the angle 9 is 

cos6> = — ^. (18) 

The coherence condition Acj) = is valid independent of Az. This is a crucial 
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property specific for Cherenkov radiation onk>Q. Thus the change of color at 
emission vertices is not important if one considers a particular a-th compo- 
nent of color fields produced at Cherenkov angle. Therefore the fields (<& a , A a ) 
and the classical current for in-medium gluodynamics can be represented by 
the product of their electrodynamical expressions (<£, A) and the color matrix 
T a . As a result, one can neglect the "rotation" of color at emission vertices 
and use in the lowest order for Cherenkov gluons the well known formulae 
for Cherenkov photons just replacing a by cxsCa for gluon currents in prob- 
abilities of their emission. Surely, there is radiation at angles different from 
the Cherenkov angle (I18p . For such gluons one should take into account the 
coherence length and color rotation considering corresponding Wilson lines 

Let us remind the explicit Abelian solution for the current with velocity 
v along z-axis 

j(r, t) = vp(r, t) = 47r 5 v5(r - vt). (19) 

In the lowest order the solutions for scalar and vector potentials are related 
so that 

AW(r,t) =ev$W(r,t), (20) 

where the superscript (1) indicates the solutions of order g. 

Therefore the explicit expressions for $ suffice. Using the Fourier trans- 
form, the lowest order solution of Eq. (llip with account of (I19p can be cast 
in the form 

•"M = < 21 > 

The integration over the angle in cylindrical coordinates gives the Bessel 
function Jo(k±r±). Integrating over the longitudinal component k z with ac- 
count of the poles due to the denominatoid and then over the transverse one 
k±, one gets the following expression for the scalar potential 



$«(r,t) = 29 ^ ~ Z ~ (22) 



yj (vt — z) 2 — r^jev 2 — 1) 



Here r± = \/ x 2 + y 2 is the cylindrical coordinate, z is the symmetry axis. 
The cone 

z = vt — r_i_v ev 2 — 1 (23) 



4 The requirement for A(j> to be a multiple of 2ir (or a weaker condition of being less or of 
the order of 1) in cases when Cherenkov condition is not satisfied imposes limits on the effective 
radiation length as it happens, e.g., for Landau-Pomeranchuk or Ter-Mikaelyan effects. 

5 These poles are at work only for Cherenkov radiation! 
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determines the position of the shock wave due to the ^-function in Eq. (|22p . 
The field is localized within this cone. The Descartes components of the 
Poynting vector are related according to Eqs. (|22p . (|20p by the formulae 



S 



s z 



(z — vt)x 



s z 




(24) 



y ~ 



so that the direction of emitted gluons is perpendicular to the cone (j23|) and 
defined by the Cherenkov angle 



which coincides with (|18|) . 

The higher order terms (g 3 ...) can be calculated using Eqs. (|15j) . (116|) . 

The expression for the intensity of the radiation is given by the Tamm- 
Frank formula (up to Casimir operators) 



It is well known that it leads to infinity for constant e. The w-dependence of 
e (its dispersion) usually solves the problem. For absorbing media e acquires 
the imaginary part. The sharp front edge of the shock wave is smoothed. The 
angular distribution of Cherenkov radiation widens. The 5-function at the 
angle (|18p . (|26p is replaced by the Breit-Wigner shape [15] with maximum 
at the same angle (but |e| in place of e) and the width proportional to the 
imaginary part. Without absorption, the potential (|22|) is infinite on the 
cone. With absorption, it is finite everywhere except the cone vertex and 
is inverse proportional to the distance from the vertex. For low absorption, 
the field on the cone increases as (Ime) _1//2 (see |29|). Absorption induces 
also longitudinal excitations (chromoplasmons) which are proportional to the 
imaginary part of e and usually small compared to transverse excitations. 
The magnetic permittivity is easily taken into account replacing e by e/i in 
the Breit-Wigner formula. 

In electrodynamics the permittivity of real substances depends on oj. More- 
over it has the imaginary part determining the absorption. E.g., Re e for water 
(see [16]) is approximately constant in the visible light region (y^e ~ 1.34), 
increases at low uj and becomes smaller than 1 at high energies tending to 1 
asymptotically. The absorption (Ime) is very small for visible light but dra- 
matically icreases in nearby regions both at low and high frequencies. Theo- 
retically this behavior is ascribed to various collective excitations in the water 
relevant to its response to radiation with different frequencies. Among them 



tan 2 (9 




= ev 




(25) 




(26) 
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the resonance excitations are quite prominent (see, e.g., [IZ])- Even in elec- 
trodynamics, the quantitative theory of this behavior is still lacking, however. 
Then, what can we say about the nuclear permittivity? 

4 The nuclear permittivity 

The partons constituting high energy hadrons or nuclei interact during the 
collision for a very short time. Nevertheless, there are experimental indica- 
tions that an intermediate state of matter (CGC, QGP, nuclear fluid ...) is 
formed and evolves. Those are J /^-suppression, jet quenching, collective flow 
(•U2), Cherenkov rings of hadrons etc. They show that there is collective re- 
sponse of the nuclear matter to color currents moving in it. Unfortunately, 
our knowledge of its internal excitation modes is very scarce, much smaller 
than in electrodynamics. 

The attempts to calculate the nuclear permittivity from first principles are 
not very convincing. It can be obtained from the polarization operator. The 
corresponding dispersion branches have been computed in the lowest order 
perturbation theory [THJ HH [20]. Then the properties of collective excita- 
tions have been studied in the framework of the thermal field theories (for 
review see, e.g., [21]). Their results with additional phenomenological ad hoc 
assumption about the role of resonances were used in a simplified model of 
scalar fields [6j to show that the nuclear permittivity can be larger than 1 that 
admits Cherenkov gluons. 

Let us stress the difference between these approaches and our considera- 
tion. In Refs. |18[ \19\ [20], [2T] the medium response to the induced current is 
analyzed. Namely it determines the nuclear permittivity. The permittivity is 
the internal property of a medium. Its quantitative description poses prob- 
lems even in QED. It becomes more difficult task in QCD where confinement 
is not understood. Therefore we did not yet attempt to compute the nuclear 
permittivity and introduced it purely phenomenologically in analogy to in- 
medium electrodynamics. Our main goal is to study the medium response to 
the external color current. Cherenkov effect is proportional to g 2 according 
to Eq. (I22p if e is constant or chosen purely phenomenologically. However 
e should tend to 1 at small g and Cherenkov effect disappears. Thus it is 
of the order of g 4 at small g. Mach waves in hydrodynamics [22] are of the 
same order. When the current j( 3 ) is treated as external one in equations of 
in- vacuum gluodynamics [U EJ [10] the effect is proportional to g G . 

We prefer to use the general formulae of the scattering theory [23] to 
estimate the nuclear permittivity. It is related to the refractive index n of the 
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medium: 

e = n 2 (27) 

and the latter one is expressed [23] through the real part of the forward scat- 
tering amplitude of refracted quant£@ReF(0°,£) 

Ren(E) = 1 + An R = 1 + ^^R e F(E) = 1 + ^^<j{E)p(E). (28) 

E A 4ftE 

Here E denotes the energy, v is the number of scatterers within a single nu- 
cleoli, 77X71- the pion mass, o{E) the cross section and p(E) the ratio of real to 
imaginary parts of the forward scattering amplitude F(E). Thus the emis- 
sion of Cherenkov gluons is possible only for processes with positive HeF(E) 
or p(E). Unfortunately, we are unable to calculate directly in QCD these 
characteristics of gluonqj and have to rely on analogies and our knowledge of 
properties of hadrons. The only experimental facts we get about this medium 
are brought by particles registered at the final stage. They have some fea- 
tures in common which (one can hope!) are also relevant for gluons as the 
carriers of the strong forces. Those are the resonant behavior of amplitudes 
at rather low energies and positive real part of the forward scattering ampli- 
tudes at very high energies for hadron-hadron and photon-hadron processes as 
measured from the interference of the Coulomb and hadronic parts of the am- 
plitudes. Re-F(0°, E) is always positive (i.e., n > 1) within the low-mass wings 
of the Breit-Wigner resonances. This shows that the necessary condition for 
Cherenkov effects n > 1 is satisfied at least within these two energy intervals. 
This fact was used to describe experimental observations at SPS, RHIC and 
cosmic ray energies. The asymmetry of the /3-meson shape at SPS [24\ and 
azimuthal correlations of in-medium jets at RHIC [U [25] were explained by 
emission of comparatively low-energy Cherenkov gluons [26\ [27] . The parton 
density and intensity of the radiation were estimated. In its turn, cosmic ray 
data [28] at energies corresponding to LHC ask for very high energy gluons to 
be emitted by the ultrarelativistic partons moving along the collision axis [5]. 
Let us note the important difference from electrodynamics where n < 1 at 
high frequencies. For QGP the high-energy condition n > 1 is a consequence 
of its instability. 

The dispersion (^-dependence of n) was taken into account. Otherwise the 
intensity of the radiation given by Eq. (126ft diverges. It can be easily incor- 
porated in Eqs. (|15p . (|16p (more precisely, in their Fourier components). The 
formula (I28h valid for n — 1 <C 1 is generalized to Lorenz-Lorentz expression 



In electrodynamics these quanta are photons. In QCD those are gluons. 
7 We can only say that KeF(E) oc g 2 at small g that confirms above estimates. 
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for larger n. The imaginary part of e can be easily accounted. In principle, it 
may be estimated from RHIC data (see |29|). 

Up to now we did not discuss one of the most important problems of the 
coordinate system in which the permittivity is defined. 

5 The rest system of the nuclear matter 

The in- medium equations are not Lorentz-invariant. There is no problem 
in macroscopic electrodynamics because the rest system of the macroscopic 
matter is well defined and its permittivity is considered there. For collisions 
of two nuclei (or hadrons) it asks for special discussion. 

Let us consider a particular parton which radiates in the nuclear matter. It 
would "feel" the surrounding medium at rest if momenta of all other partons 
(or constituents of the matter), with which this parton can interact, sum to 
zero. In RHIC experiments the triggers which registered the jets (created by 
partons) were positioned at 90° to the collision axis. Such partons should be 
produced by two initial forward-backward moving partons scattered at 90°. 
The total momentum of other partons (medium spectators) is balanced be- 
cause for such geometry the partons from both nuclei play a role of spectators 
forming the medium. Thus the center of mass system is the proper one to con- 
sider the nuclear matter at rest in this experiment. The permittivity must be 
defined there. The Cherenkov rings consisting of hadrons have been registered 
around the away-side jet which traversed the nuclear medium. This geometry 
requires however high statistics because the rare process of scattering at 90° 
has been chosen. 

The forward (backward) moving partons are much more numerous and 
have higher energies. However, one can not treat the radiation of such a pri- 
mary parton in c.m.s. in the similar way because the momentum of spectators 
is different from zero i.e. the matter is not at rest. Now the spectators (the 
medium) are formed from the partons of another nucleus only. Then the rest 
system of the medium coincides with the rest system of that nucleus and the 
permittivity should refer to this system. The Cherenkov radiation of such 
highly energetic partons must be considered there. That is what was done 
for interpretation of the cosmic ray event in [5] . This discussion clearly shows 
that one must carefully define the rest system for other geometries of the 
experiment with triggers positioned at different angles. 

Thus our conclusion is that the definition of e depends on the experiment 
geometry. Its corollary is that partons moving in different directions with 
different energies can "feel" different states of matter in the same collision 
of two nuclei because of the dispersive dependence of the permittivity. The 
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transversely scattered partons with comparatively low energies can analyze the 
matter with rather large permittivity corresponding to the resonance region 
while the forward moving partons with high energies would "observe" low 
permittivity in the same collision. This peculiar feature can help scan the 
(lnx, Q 2 )-plane as it is discussed in |30j. It explains also the different values 
of e needed for description of RHIC and cosmic ray data. 

These conclusions can be checked at LHC because both RHIC and cosmic 
ray geometry will become available there. The energy of the forward moving 
partons would exceed the thresholds above which n > 1. Then both types 
of experiments can be done, i.e. the 90°-trigger and non-trigger forward- 
backward partons experiments. The predicted results for 90°-trigger geometry 
are similar to those at RHIC. The non-trigger Cherenkov gluons should be 
emitted within the rings at polar angles of tens degrees in c.m.s. at LHC by 
the forward moving partons (and symmetrically by the backward ones). This 
idea is supported by some events observed in cosmic rays [28j |2T]. 

6 Conclusions 

The equations of in- medium gluodynamics (|15[) . (|16p are proposed. They re- 
mind the in-medium Maxwell equations with non-Abelian terms added. Their 
lowest order classical solutions are similar (up to the trivial color factors) 
to those of electrodynamics (I22p . especially, for Cherenkov gluons. The nu- 
clear permittivity of the hadronic medium is related to the forward scattering 
hadronic amplitudes and its possible generalization is discussed. This def- 
inition asks for the distinction between the different coordinate systems in 
which the Cherenkov radiation (and nuclear permittivity) should be treated 
for partons moving in different directions with different energies. 

This consideration has led to explanation of several effects observed at 
SPS, RHIC, cosmic ray energies and predicts new features at LHC |30j. Some 
estimates of properties of the nuclear matter formed in ultrarelativistic heavy- 
ion collisions have been done and are predicted. 
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